Characterizing topological order by studying the ground states of an infinite cylinder 
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Given a microscopic lattice Hamiltonian for a topologically ordered phase, we describe a tensor 
network approach to characterize its emergent anyon model and, in a chiral phase, also its gapless 
edge theory. First, a tensor network representation of a complete, orthonormal set of ground states 
on a cylinder of infinite length and finite width is obtained through numerical optimization. Each 
of these ground states is argued to have a different anyonic flux threading through the cylinder. In 
a chiral phase, the entanglement spectrum of each ground state is seen to reveal a different sector 
of the corresponding gapless edge theory. A quasi-orthogonal basis on the torus is then produced 
by chopping off and reconnecting the tensor network representation on the cylinder. Elaborating 
on the recent proposal of [Y. Zhang et al. Phys. Rev. B 85, 235151 (2012)], a rotation on the 
torus yields an alternative basis of ground states and, through the computation of overlaps between 
bases, the modular matrices S and U (containing the mutual and self statistics of the different anyon 
species) are extracted. As an application, we study the hard-core boson Haldane model by using the 
two-dimensional density matrix renormalization group. A thorough characterization of the universal 
properties of this lattice model, both in the bulk and at the edge, unambiguously shows that its 
ground space realizes the v — 1/2 bosonic Laughlin state. 



For many decades, characterizing the emergent order of 
an interacting quantum many-body system from its mi- 
croscopic description has been regarded as an extremely 
challenging task. Consider, for instance, a lattice Hamil- 
tonian H suspected of realizing some form of topological 
order H , say a given Laughlin state Q or a quantum spin 
liquid 0, IH — states of considerable interest both in the 
study of exotic many-body phenomena such as the frac- 
tional quantum Hall (FQH) effect |5j or high-temperature 
superconductivity 0t3? an ^ m ^ ne design of a quantum 
computer based on topological protection to decoherence 
[9|, [Toj| . Due to a lack of theoretical and computational 
tools, assessing whether the low energy sector of H is 
indeed topologically ordered, and then establishing what 
type of topological order it realizes, have traditionally 
been considered very difficult problems. However, recent 
advances in the understanding of many-body entangle- 
ment have progressively brought us closer to being able 
to tackle them. On the theoretical side, new ways of di- 
agnosing the presence and type of topological order from 
knowledge of the ground state wave- function alone, based 
on entanglement entropy pjKHj, entanglement spectrum 
[l3| and modular transformations [3, EH , have been put 
forward. On the computational side, the advent of ten- 
sor networks [3^13 makes it now possible, by mimicking 
the structure of entanglement, to efficiently represent a 
large class of low energy many-body states. 

In this paper, by combining and building upon the 
above developments, we describe an approach for char- 
acterizing the topological order emerging from a micro- 
scopic Hamiltonian H in a two dimensional lattice. The 
key of our approach is the computation of a basis of 
ground states of i7, first on an infinite cylinder, and then 
on a finite torus. Each ground state has a flux through 
the cylinder /torus corresponding to one of the charges 
of the emergent anyon model. These ground states are 
encoded in a tensor network representation, from which 



we can extract the universal properties of the emergent 
edge and bulk theories. For concreteness, we focus on a 
sjDecinc lattice Hamiltonian, namely the Haldane model 
(2^ on the honeycomb lattice with hard-core bosons [25| , 
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where we set (j> = 0.4tt and (*,£',£") = (1,0.6,-0.58), 
where t, t' and t" stand for the hoping amplitudes be- 
tween nearest, next-nearest and next-next-nearest neigh- 
bors, respectively. For these parameters, Ref. [25| found 
two quasi-degenerate ground states and a non-trivial 
Chern number C = 1 using exact diagonalization on 
small tori, while Ref. [22| computed a non-zero topolog- 
ical entanglement entropy (TEE) [O, 03 7 ~ log y/2 by 
studying finite cylinders with the density matrix renor- 
malization group (DMRG) [16|. These results are sug- 
gestive of a v = 1/2 bosonic FQH state [2|, which has a 
chiral semion [26| in the bulk and a SU (2)i Wess-Zumino- 
Witten (WZW) conformal field theory (CFT) [27| at the 
edge. Here we will not make use of any previous knowl- 
edge about i^Hai., and yet will produce a detailed char- 
acterization of both bulk and edge theories. Our ap- 
proach can be readily applied to other Hamiltonians on 
honeycomb, triangular and kagome lattices (with tt/3 ro- 
tational symmetry), and can be generalized to arbitrary 
lattices (even without any rotational symmetry) [28[ . 

Complete set of ground states on an infinite cylinder. — 
Consider a local Hamiltonian H on a cylinder of size L x x 
L yi with infinite length L x = 00 and finite width L yi see 
Fig. [Ua). If the width L y was also infinite, H would have 
a finite number N of exactly degenerate ground states, 
one for each species i of anyon in the emergent anyon 
model [lj. A finite L y implies, however, that the ground 
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FIG. 1. (Color online) (a) Honeycomb lattice on a cylinder 
with L x = oo and L y = 4 (as measured in number of unit 
cells). A blue, dashed line indicates the entanglement cut 
used in dividing the cylinder into two semi-infinite cylinders, 
(b) Ground states |^ yl ) and |^ yl ) of H Ua i. in Eq. [Qwith flux 
i = 1,8. (c) Corresponding density matrices pi and p s for the 
left half of the cylinder. 

state degeneracy is broken. This is due to the creation 
of virtual pairs of anyons that annihilate after describing 
a non-trivial loop in the y direction [9]. If the width 
L y of the cylinder is much larger than the correlation 
length £ in the ground states of H, L y ^> £, we expect 
the topological order to be still (approximately) realized, 
and the "low energy" space of H to still decompose into 
TV distinct topological sectors. Each of these sectors is 
characterized by a ground state |^- y )' 

H\9?)=E i \9? 1 ), (2) 

which, crucially for the present approach, has an anyon 
of type i propagating inside the cylinder in the x direc- 
tion, see Fig. QJb) - that is, |\fr- yl ) is an eigenvector of 
the dressed Wilson loop operator [32j encircling the cylin- 
der in the y direction. This is so because virtual pairs 
of anyons encircling the cylinder in the x direction only 
renormalize the energy E{ of |\I/J yl ) (see also Appendix A). 
[It can be further argued, using degenerate perturbation 
theory, that the gap A^- = E{— Ej between ground state 
energies is exponentially suppressed with L y but exten- 
sive in L x , Aij « e~ aLy L x , and thus infinite. An infinite 
gap does not jeopardize the present approach.] 

Tensor networks. — Each ground state |^J yl ) can be 
represented using a tensor network with a finite unit cell 
of tensors that is repeated throughout the infinite cylin- 
der (see Appendix B). A natural choice is to use a pro- 
jected entangled-pair state (PEPS) [ItJ , optimized e.g. 
through energy minimization [I?} or imaginary time evo- 
lution [19|. The incurred computational cost is at most 
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FIG. 2. (Color online) (a) MPS for a state |^ yl ) of an in- 
finite cylinder with L y = 4, see Fig. [TJa). The bond in- 
dices that connect tensors form a snake. The MPS unit cell 
is made of eight tensors that are repeated throughout the 
tensor network, (b) MPS for a state \^l° T ) °f a torus with 
L x x L y x 2 = 32 sites, see Fig. SJa), obtained by reconnect- 
ing four MPS unit cells of the state |^- yl ). 



linear in L y , which allows for the constraint L y ^> £ to 
be easily met. In the present application we study the 
ground states of H Hald _ using instead a matrix product 
state (MPS) optimized using DMRG 0, see Fig. [2{a). 
The computational cost grows now exponentially with 
the width L y of the cylinder, and therefore only small 
values of L y can be afforded. However, if the condition 
L y ^> £ can still be met, DMRG offers a more accurate 
and dependable implementation of our approach. 

Specifically, with an infinite DMR G_|29l l30| algorithm 
in two dimensions (detailed in Ref. |28j), we considered 
infinite cylinders with width L y = 4,6 and 8 (as mea- 
sured in number of unit cells, see Fig. [1]). For each 
value of L y , we repeated the optimization of an infinite 
MPS several times. Each optimization produced one of 
two states, with each state occurring roughly half of the 
runs. We thus conclude that there are at least two ground 
states or, equivalently, that the emergent anyon model 
has N > 2 anyon charges. Let us denote these two ground 
states |\E^ yl ) and |^s yl ), anticipating that they correspond 
to anyon types i = 1, or identity charge, and i = s, or 
semion charge, respectively, to be confirmed later. A 
very small correlation length £i « £ s « 0.52 for both 
states verifies that L y ^> £ and thus we have effectively 
reached the thermodynamic limit. 

Quantum dimensions from entanglement entropy. — 
The tensor network representation of state |\E^ yl ) readily 
provides access to the eigenvalues {pi, a } of the reduced 
density matrix pi for half of the infinite cylinder, see 
Fig. HJc). From the scaling of the entanglement entropy 
S(pi) = — ^2 a Pia log(pia) as a function of the width L y , 
namely Si(L y ) = aL y — ji, we can estimate the TEE 7^ 
(Hill for ground state |^ yl ). We obtain 71 w 0.3455, 
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FIG. 3. (Color online) Entanglement spectrum of the reduced 
density matrix pi for half an infinite cylinder for ground state 
l^ 1 ), for i = 1 and s, Fig. Ufa). The vertical axis shows 
Ei i(X = — log(pi, a ) (up to global shift and re-scaling), where 
{pi,a} are the eigenvalues of pi. The horizontal axis shows 
the momentum of the corresponding eigenvector of pi, which 
we artificially extend beyond its 2tt periodicity. A tilted line 
connects all eigenvalues of pi with the same particle number 
M, also indicated (notice that ifnai. has a U(l) symmetry 
corresponding to particle number conservation). We find a 
degeneracy pattern {1,1,2,3,5,---} (according to the parti- 
tion numbers) along each tilted line of fixed M, characteristic 
of a bosonic Gaussian theory. 



sponds to the chiral boson vertex operator e 1 ^!^ 2 and its 
descendants, which justifies our previous identification 
of 1^7 ) an d |^s yl ) with the anyonic charges i = 1 and 
s, respectively. Thus, from the states {|\I/J y1 )} we have 
unambiguously identified the edge CFT of the emergent 
anyon model. 

From an infinite cylinder to a finite torus. — Using an 
infinite cylinder has the major advantage, compared to 
previous DMRG studies on a finite cylinder |2ll-[23| , that 
it provides a complete set of ground states {|^- yl )i °f a 
local Hamiltonian H. In addition, we can also produce a 
complete basis {|\I/^ or )} for the (quasi-degenerate) ground 
space of H on a finite torus of size L x x L yi where the 
choice L x = L y ensures that also L x ^> £. This is ac- 
complished by reconnecting a region of size L x x L y of 
the tensor network for {|^ y1 )} into a torus, see Fig. [2fb) 
(and Appendix B). Notice that state |^* or ) has an anyon 
of type i inside the torus. Following Ref. [14[ we can 
use the basis j|\E^ or )i to characterize the emergent anyon 
theory of H [331 ] . 

Bulk anyon model from modular transformations. — 
Our goal is to compute the modular matrices S and U 
characterizing the mutual and self statistics of the emer- 
gent anyon model, 



7 S ~ 0.986 x 7i, very similar to the finite cylinder result 
of Ref. [22|] for one ground state. The expression [ll|, [3lj 



■log(di/D), 



D 



(3) 



allows us to compute each quantum dimension di and the 
total quantum dimension D (see Appendix D). First, we 
note that the sum ^^di/ D) 2 only reaches unity when 
it includes all the anyon charges. We can therefore use 
this sum to establish whether we have obtained a com- 
plete set of ground states. Our estimates of 7^ yield 
(di/D) 2 + (d s /D) 2 w 1.007. Thus, i7 Hal . only has two 
ground states, i.e. the emergent anyon model has N = 2 
types of anyons [25|. Then, since the identity charge al- 
ways exists and has quantum dimension d\ = 1, from 
7 S — 71 = —\og(d s / d\) we find d s = 1.005, and from 
71 = — \og(d\/D) we find D = 1.413, very close to the 
exact d s = 1 and D = y/2 « 1.4142 of the semion model 

Edge chiral CFT from entanglement spectrum. — The 
entanglement spectrum (ES) [13[ for |^- yl ) (see Appendix 
C), plotted in Fig. [3j reveals that each ground state of 
^Hai. is associated with one of the two primary field of 
the chiral 517(2) 1 WZW CFT, namely the identity and 
chiral boson vertex operator e 1 ^!^ 2 (see e.g. Chap 15.6 
of Ref. [27j), which are seen to be a singlet and a dou- 
blet of an emergent SU(2) symmetry. Specifically, the 
ES of \^l yl ) is organized according to the scaling dimen- 
sions and conformal spins of the identity primary field 
of this CFT and the tower of all of its (Kac-Moody and 
Virasoro) descendants, whereas the ES of \^ c s yl ) corre- 
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Here, entry Sij determines the phase acquired by an any- 
onic charge i when encircling an anyonic charge j; c is the 
central charge of the anyon model; and 6i corresponds to 
the phase acquired when an anyonic charge i is exchanged 
with another identical anyonic charge i (and thus 0\ = 1). 
From Su = Su = di/D and Y^i(Su) 2 = 1 (5 is a uni- 
tary matrix) we can also extract the quantum dimensions 
di and total quantum dimension D, whereas the fusion 
rules i x j — Nij k k follow from the Verlinde formula 
Nij k = ^2 rn (Si m Sj rn (S rn ky I Si m ) [34J- 

A 7r/3 rotation R^/3 on the torus, see Fig. [H cor- 
responds to applying the modular transformation us -1 , 
where u and 5 generate the group of modular transfor- 
mations on a torus with defining vectors w\ and 1^2, 



s : 



Following Ref. Q, the overlaps Vij = (^T r \ R n/s\^T^ 
between the bases {|^- or )l and {^7r/3l^j° r )} form tne 
unitary transformation V = D^US~ 1 D, where D is a 
diagonal matrix of phases Djj = corresponding to a 
phase freedom in choosing |^ or ), and S and U generate a 
representation of the modular group on the ground space 
of i7 H ai - We thus have 
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FIG. 4. (Color online) (a) Torus with L x — L y — 4, made 
of L x x L y x 2 = 32 sites, and with generating vectors u>i 
and W2- Four sites are marked as o,p, r for reference, (b) 
Torus obtained after a ty/3 rotation R n /3, corresponding to a 
modular transformation that maps (wi,W2) into (w^w'2) = 
(wi —W2,wi). (c) By using the periodicity by vectors (wi,W2) 
[equivalent ly, (w[,w 2 )], the sites in torus (b) can be mapped 
back into torus (a). This allows us to compare states defined 
on the two tori and compute V in Eq. \E\ following Ref. [IH ]. 



where in the second line we use that Su,Sn > and 
61 = 1. Building upon Ref. 0, we note that from ma- 
trix V we can actually compute both U and S. Indeed, 
from Vn we obtain the central charge c and Su; then 
from Vi B we obtain e l ^ s ~^^ and S s \\ then from V a \ we 
obtain 6 8 and 5i s ; finally, V 88 yields S ss . [Similarly, for an 
arbitrary anyon model, from the scalar products Vij we 
can completely determine both Sij and Uij , see Appendix 
E] . The exact evaluation of Vij from a tensor network rep- 
resentation turns out to have a very large computational 
cost. We compute Vij approximately using Monte Carlo 



sampling on a periodic MPS [35|. For L y = 6 (torus with 
6 x 6 x 2 = 72 sites), we obtain 
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with sampling noise on the order of 10 3 , which accu- 
rately reproduce the exact S and U matrices of a chiral 
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semion model, namely 

In particular, from Eq. [6] we can accurately extract the 
quantum dimensions d\ = d s = 1 and total quantum di- 
mension D = y/2 (independently of the previous entropy 
calculations) and the Z2 fusion rules 1 x 1 = s xs = 1 
and 1 x s = s x 1 = s; whereas from Eq. [7] we determine 
a central charge c = 1 (modulus 24) and the twist S = i 
characteristic of a semion. 

In this paper we have shown how to produce, start- 
ing from a microscopic lattice Hamiltonian H, a detailed 
characterization of the emergent topological order. A 
key step is to obtain a tensor network representation for 
a complete set of ground states of H, first on an infi- 
nite cylinder and then on a finite torus. Here we have 
considered several ground state properties that can be 
extracted directly from the optimized tensor network. 
With further manipulation, it is also possible to obtain 
an explicit representation of quasi-particle excitations of 
H, including those with fractional quantum numbers, as 
discussed elsewhere [28j. 

The authors thank Jaume Gomis and Xiao- Gang Wen 
for guidance in identifying the edge theory resulting 
from the numerical simulations, and Oliver Buerschaper, 
Tarun Grover, Todadri Senthil, and Ashvin Vishwanath 
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Appendix A: Ground states on an infinite cylinder 
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Let H be a local lattice Hamiltonian corresponding 
to a topologically ordered phase, and let £ be the largest 
correlation length in the ground space of H on an infinite 
lattice. In this appendix we argue that, when placed on 
a torus of size L x x L y with L x = oo and L y ^> £, each of 
the (quasi-degenerate) ground states of H has a 

well-defined anyon flux i in the x direction, as measured 
by a dressed [32] Wilson loop operator T y encircling the 
torus in the y direction, see Fig. [5l Then we argue that 
when placed on a cylinder of size L x x L y with L x = oo 
and L y ^> £ (and convenient boundary conditions to be 
specified below), H also has analogous (quasi-degenerate) 
ground states {|\I/^ y1 )}. This last set of ground states 
can be represented and computed with tensor network 
techniques, as discussed in the main text of this paper 
and in Appendix B. 



1. Ground states of a fixed-point Hamiltonian 

Let us first consider a fixed-point Hamiltonian H ^ de- 
scribing the same phase as H (but with vanishing corre- 
lation length, £ = 0). We will assume that V = H-HW 
is a small perturbation that can be expressed as a sum of 
local terms and treated using (degenerate) perturbation 
theory. 

Let denote a basis of ground states of 

on a torus of size L x x L y , where state l^f^) has anyon 
flux i in the x direction, i.e., I^f^) is an eigenvector 
of the Wilson loop operator T y ^ encircling the torus in 
the y direction, see Fig. [5fa)-(b). Clearly, the states 
{|^°^)} are just one possible basis in the ground space 
of H( \ but this particular choice of basis plays a dis- 
tinguished role below. We remind the reader that all 
ground states of the topologically ordered Hamiltonian 
are locally indistinguishable - their reduced density 
matrices are identical on any contractible region - and 
that in order to transform the ground state |^-°^) into 

another ground state |^°^), a non-local operator closing 
a non-trivial loop in the x direction, such as a Wilson 
loop operator r^°\ is required. 

Let us now assume that, starting with the lattice in 
the ground state l^- ^) of H^°\ we adiabatically switch 
on perturbation V, so that Hamiltonian is changed from 
to H. Let \S&i) denote the resulting state of the 
system, 
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FIG. 5. (Color online) (a) Torus of size L x x L y . A Wilson 
loop operator Ty encercling the torus in the y direction mea- 
sures the anyon flux in the x direction, (b) The same torus, 
with anyonic flux i in the x direction, (c) Modifying a fixed- 
point Hamiltonian with a perturbation V has the effect, 
in perturbation theory, of creating virtual pairs of anyons that 
propagate through the torus. The figure represents a pair of 
anyons j and j that annihilate after closing a non-trivial loop 
in the y direction. This may shift the energy of \^i) by an 
amount that depends on the anyon type i. (d) For L x = oo 
and finite L y , the ground state degeneracy is broken in the 
\^i) basis. 



For a generic choice of sizes L x and L y , will no longer 
be an eigenvector of H. However, we will argue below 
that for L x = oo and finite L y ^> £, |\f^) is still an 
eigenvector of H. In preparation, we first consider a torus 
with (L x ,L y ) = (oo, oo). 



2. Ground states of if on a torus with L x = L y = oo 

For L x = L y = oo, let us analyze the effect of pertur- 
bation V (which by hypothesis does not break the topo- 
logical order) using perturbation theory. Perturbation V 
can be interpreted as to allow for the creation of virtual 
pairs of anyons that propagate throughout the torus @ . 
If one virtual pair annihilates after closing a non-trivial 
loop, this may change the topological properties of the 
ground state. However, we notice that, due to the infi- 
nite size of the torus in any direction, no virtual pair of 
anyons can close a non-trivial loop at m th order in pertur- 
bation theory, for any finite m, because the m th power of 
V decomposes as the sum of terms each of which still has 
support on a contractible region only. In other words, 
for any finite m, the overlap (^f^ |(V) m |^^) is equal to 
some value ((^) m ) independent of the flux i (since all 
ground states of are identical on any contractible 



region), whereas for i ^ j the overlap (^°^|(V) m |^ u; ) 
(Al) vanishes (because (V) m cannot close a non-trivial loop 



l^ (0) ) 



l*i>- 



in the x direction), 

<*i 0) |(y) m |*f } ) = Sij((V) m ). (A2) 

As a result, in degenerate perturbation theory each 
state \^i) is obtained from the corresponding l^- ^) by 
means of the same local modifications, so that all the 
states \S&i) remain locally indistinguishable, and form a 
complete set of (exactly degenerate) ground states of H 
with renormalized energy E, i.e. H\$i) = E\^i), and 
finite correlation length £. In particular, state still 
has anyon flux i in the x direction, as measured now by 
a dressed [32| Wilson loop operator T y in the y direction. 

3. Ground states of H on a torus with L x — oo and 
finite L y ^> £ 

For L x = oo and finite L y ^> £, the ground state de- 
generacy of H is broken in perturbation theory due to the 
fact that virtual pairs of anyons can now annihilate with 
each other after their combined trajectories have closed a 
non-trivial loop in the y direction, see Fig. [5jc), produc- 
ing an energy contribution ((V) m )i = (^ ( ° } |(V) m |^ (0) ) 
(occurring first for some finite m proportional to L y ) that 
may now depend on the flux i present in the x direction 
of the cylinder. On the other hand, no virtual loop can 
still be closed in the x direction for any finite m, since 
L x = oo, and thus ) |(V r ) m |*5° ) > = for j, or 

(M 0) \(vr\*? ) ) = kj((vr)i- (A3) 

As a result, the ground state degeneracy is broken pre- 
cisely in the basis, so that each |\E^) remains an 
exact eigenvector of H with (now distinct) energy E{, 
H\^i) = Ei\Vi). A gap Aij = Ei- Ej opens. Since this 
first occurs in perturbation theory at a finite order m that 
is proportional to L y , the gap decays exponentially in L y , 
w e ~ aL v 5 for some constant a [9]. However, tunnel- 
ing of virtual anyon pairs in the y direction can occur 
anywhere in the x direction, and therefore the gap is also 
proportional to L x , « e~ aLy L x , and thus infinite. 
[This will turn out to be of no relevance in our numerical 
simualtions based on energy minimization, because the 
simulation still gets trapped in some topological sector, 
independently of the gaps A^]. Importantly, once again 
the state |\f^) has anyon flux i in the x direction, as mea- 
sured now by a dressed [32| Wilson loop operator T y in 
the y direction. We continue to refer to state |\E^) as a 
ground state of H, since it is its ground state in the i flux 
sector. 



4. Local equivalence between the ground states of 
H on a torus and a cylinder when L x — oo and finite 

Ly >£• 

We have argued above that the combined effect of per- 
turbation V and finite L y <C £ on the torus is that Hamil- 
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FIG. 6. (Color online) (a) Schematic representation of a 
ground state |^) of H on a torus of size L x x L y with 
L x — oo and finite L y , with anyon flux i in the x direction, (b) 
Schematic representation of an analogous ground state |^^ yl ) 
of H on a cylinder of the same size L x x L y , where the anyon 
flux is created and annihilated at the x — =boo boundaries. 



tonian H has eigenstates (referred to as ground 

states) with well defined anyon flux i in the x direction. 

Let us now consider a cylinder of size L x x L y , with 
= oo and finite L y ^> £, and a Hamiltonian on this 
cylinder that decomposes into: (i) a bulk Hamiltonian 
that is locally identical to the Hamiltonian H above for 
the torus, and (ii) boundary Hamiltonian terms H~°° 
and H + °° for the x = — oo and x = oo boundaries of the 
cylinder. [Here L x = oo describes the limit of a cylinder 
with finite L x , for which an analogous Hamiltonian can 
be defined]. We assume that H ±oc is such that an anyon 
of type i can be created at one boundary and annihilated 
at the other boundary at a finite energy cost. Then the 
ground states of H cy \ denoted {|^- y1 )}, have anyon flux 
i through the cylinder, see Fig. [6j as can be argued 
in a similar way as we discussed above for a torus. In 
addition, on any finite section of the cylinder |\E^ yl ) is 
identical to the analogous ground state |\E^) on the torus, 
since these states only have short correlations and result 
from minimizing a local Hamiltonian, which is the same 
on the torus than in the bulk of the cylinder. 

The particular choice of the above boundary Hamilto- 
nian terms H ±OQ is actually of no practical importance 
for our numerical approach. Indeed, in a tensor network 
computation of the ground states of an infinite cylinder, 
as used in this paper, the optimization of the bulk tensors 
will not take into account energy contributions coming 
from the boundary. 



Appendix B: Tensor network representations on an 
infinite cylinder 

In this work we have used a tensor network to represent 
states of a cylinder of size L x xL y , where L x = oo and L y 
is finite. We consider two different tensor networks: the 
projected entangled-pair state (PEPS) and the matrix 
product states (MPS). 
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FIG. 7. (Color online) (a) PEPS representation for a state 
|^ cyl ) on a cylinder of size (L x , L y ) = (oo,4) (honeycomb lat- 
tice). Notice that the PEPS consists in repeating two tensors 
A and B, and thus the PEPS unit cell coincides with the 
natural unit cell on the honeycomb lattice, (b) MPS repre- 
sentation for the same state. Here the MPS unit cell is made 
of eight different tensors {A, B, • • • , H}, connected in a one- 
dimensional geometry that snakes around the cylinder. Each 
MPS unit cell corresponds to four (in general, L y ) unit cells 
of the honeycomb lattice. 



1. PEPS on an infinite cylinder 



der, which allows us to consider large values of L 



2. MPS on an infinite cylinder 

An alternative tensor network representation for states 
of an infinite cylinder is a matrix product state (MPS), 
which is the basis of the density matrix renormalization 
group (DMRG) [Hj]. In the example of Fig. [7(b), a unit 
cell of tensors is also repeated throughout the cylinder, 
but the MPS unit cell is L y times the size of the natural 
unit cell on the honeycomb lattice. An MPS tensor A 
has components (A) a p s , where the bond indices a and 
(3 take \ values and the physical index s takes d values. 
Thus, the number of variational parameters in an MPS 
for an infinite cylinder is of the order of 2L y dy 2 . Key to 
this discussion is the fact that the value of \ necessary to 
obtain an accurate description of a ground state |\E^ y ) * s 
expected to grow exponentially with L y , making the MPS 
representation (and its manipulation, with cost scaling 
as x 3 ) inefficient. As a result, only systems with small 
values of L y can be afforded computationally. For states 
with a sufficiently small correlation length £ (as given, 
for instance, by the transfer matrix correlation length 
£ TM , see next subsection), such small values of L y might 
already be sufficient in order to make finite-size effects 
negligible and thus obtain the thermodynamic properties 
of the system. This is precisely the case with the ground 
states of Hamiltonian H Hal in Eq. [TJ 

Twenty years of successful DMRG experience makes 
this approach much better understood than PEPS, which 
is only a few years old. DMRG offers a more accurate and 
dependable implementation of our approach whenever it 
can be applied - that is, when the condition L y ^> £ can 
be met for sufficiently small cylinder width L y . 



A natural option is to use a projected entangled pair 
state (PEPS) [17] with periodic boundary conditions in 
the y directions. In the presence of translation invariance, 
the state of an infinite system can be represented by re- 
peating a unit cell of tensors throughout the lattice [HI . 
For instance, in a honeycomb lattice, two tensors A and 
B can be used to form a PEPS unit cell, see Fig. [TJa). In 
this example, the PEPS unit cell coincides with the nat- 
ural unit cell on the honeycomb lattice. Notice that this 
representation is very compact. Tensors A and B have 
components (A) a p^ s and (B) a p ls , where the bond in- 
dices a, /?, and 7 take D values each (e.g. a = 1, • • • , D) 
and the physical index takes d values, where d is the di- 
mension of the Hilbert space of one lattice site (d = 2 in 
Hamiltonian H Hal , since on each site we may have zero or 
one bosons). This amounts to about 2dD 3 complex co- 
efficients, or variational parameters, to characterize the 
PEPS. For a gapped phase, the value of D necessary to 
obtain an accurate description of a ground state |^^ yl ) 
is not expected to depend on L y . Importantly, the com- 
putational cost of optimizing these tensors can also be 
made essentially independent of the size L y of the cylin- 



3. Transfer matrix 

In the tensor network representation of a state of an in- 
finite cylinder, correlations in the x direction can be seen 
to be mediated by powers of a transfer matrix (TM). This 
TM is expressed as a tensor network, see Fig. [5J made of 
reduced tensors. For instance, for a PEPS on a honey- 
comb lattice with tensors A and £?, the reduced tensors 
a and b are built from A and A* (complex conjugate), 
and B and 5*, respectively, according to 

(a)(aa')(00')(77') = E^ a ^^^'7'»' ( B1 ) 
s 

where (a, a 7 ) denotes a double index. 

Let {Ai, A2, • • • } denote the eigenvalues of the TM, or- 
dered in deceasing absolute value, |Ai| > | A2 1 > • • • • In 
a normalized tensor network, the largest eigenvalue Ai 
of the TM is always one, Ai = 1, and the decay of cor- 
relators c(x) in the x direction is governed by the rest 
of eigenvalues {A2, A3, • • • }. In particular, the dominant 







FIG. 8. (Color online) (a) The PEPS transfer matrix on the 
x direction is obtained by building the reduced tensors a and 
b, see Eq. IB11 and connecting these tensors periodically in 
the y direction. The inset shows tensors a and b in terms of 
tensors A, A* , and B, B* , respectively, (b) The MPS transfer 
matrix is built in an analogous way. The inset shows two of 
the reduced tensors, a and b. The second largest eigenvalue 
A2 of the TM determines the largest correlation length of the 
system in the x direction. 




FIG. 9. (Color online) (a) PEPS representation for a state 
|^ tor ) on a finite torus of size (L x , L y ) = (4, 4) obtained by re- 
connecting the PEPS tensors A and B that form a PEPS unit 
cell for a state |^ cyl ) on a cylinder of size (L x ,L y ) = (00, 4). 
Notice that the resulting tensor network has the topology of 
a torus, (b) Similar construction for a MPS. 



correlations decay exponentially as 

c(x) « (A 2 f = e -^ TN , (B2) 

where £ TM = — l/log(A2) defines the TM correlation 
length. [Here we have assumed, for simplicity, that 
I A2 1 > I A3 1; for A2 = A3 a similar expression holds]. 

When studying topological order with infinite cylin- 
ders of finite L yi the TM correlation length £ TM of the 
ground states offers a natural consistency check for the 
assumption that the value of L y is large enough to be 
representative of the thermodynamic limit. If the size 
L y of the cylinder is much larger than £ TM , then we in- 
deed expect finite size effects to be very small, justifying 
a posteriori the use of a cylinder with that size L y . If, 
on the contrary, L y is somewhat comparable to £ TM (or 
£tm depends significantly on the size L y of the cylinder), 
it is likely that finite size effects are still dominant. 



4. From an infinite cylinder to a finite torus 

Given a tensor network representation of state |\I/ cyl ) on 
an infinite cylinder with finite size L y in the y direction, 
it is possible to obtain a tensor network representation of 
a state |\I/ tor ) on a torus with the same size L y and some 
finite size L x in the x direction by conveniently recon- 
necting into a torus the tensors in a proper subregion of 



the cylinder, as illustrated in Fig. [9] If in the cylinder we 
had £ TM <C L y for a ground state |\E^ yl ) of H with anyon 
flux i in the x direction, then we expect that the resulting 
state |\E^ or ) on a torus of size (L x ,L y ) with L x > L y be 
a good approximation to the ground state of H on that 
torus with the same anyon flux i in the x direction. The 
results of the present paper offer a clear confirmation of 
this expectation. 



Appendix C: Entanglement spectrum of half an 
infinite cylinder from a matrix product state 

Given an infinite MPS for a state |\I/ cyl ) of an infi- 
nite cylinder, see Fig. EJb), we obtain a tensor network 
for the reduced density matrix p for half of the cylin- 
der by cutting the MPS into two halves, and then trac- 
ing over the open bond index, see Fig. [TOt a)-(b). [Here 
we assume that the open index is normalized as to con- 
tain the Schmidt coefficients of that cut]. Translations 
by one lattice unit cell (or two sites) in the y direction 
generate a finite Abelian group Z^ y of transformations 
on the infinite cylinder. If is invariant under such 
translations, then so is p. Therefore each eigenvector 

of p transforms according to some phase e L y , where 
K = 0, 1, • • • , L y — 1 is the momentum of that eigen- 
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(a) 



(b) 



(c) 




FIG. 10. (Color online) (a)-(b) From an MPS for the state 
|\I/ cyl ) of an infinite cylinder, see Fig. a MPS for the re- 
duced density matrix p for half the infinite cylinder is obtain 
by tracing out the open index in this figure (assuming the 
proper normalization of the MPS tensors), (c) Alternatively, 
a density matrix p with the same spectral properties as p is 
obtained from this figure by tracing out all the physical in- 
dices (again, assuming the proper normalization of the MPS 
tensors) . 



(a) 




FIG. 11. (Color online) (a) A translation by one lattice unit 
cell (made of two sites) in the y direction is implemented by 
permuting the indices of the MPS unit cell (made of 8 ten- 
sors). If state |^ cyl ) is invariant under such a translation, 
then the MPS unit cell before and after the corresponding 
permutation of sites still describes the same state, and it is 
therfore related by a transformation R that, under proper nor- 
malization of the MPS tensors, is unitary and applies a phase 

e L y to each eigenvector of p, where K = 0, 1, • • • , L y — 1 
is the momentum of that eigenstate. (b) Matrix R is obtained 
as the fixed point of a (mixed) transfer matrix. 



vector. Our goal is to obtain the spectrum {p a ,K a } of 
P- 

If instead of tracing the open bond index, we trace out 
the infinitely many physical indices of half of the cylinder 
[now assuming that the open index does not contain the 
Schmidt coefficients of that cut] , then we obtain a density 
matrix p, living on the vector space of the open bond 
index, with the same spectral properties as p, see Fig. 
[TOl c). We can then obtain the eigenvalues {p a } of p 
by diagonalizing p. In addition, we can also obtain the 
momentum K a associated to each eigenvalue p a by the 
construction in Fig. [TTJ where the matrix R is a unitary 
representation of the Abelian group of translations 
in the y direction. 



Appendix D: Computation of quantum dimensions 
from Entanglement Entropy 

In this appendix we explain how to extract the quan- 
tum dimensions {di} and total quantum dimension D = 
V^i^) 2 ' f rom the ground states {|\I/^ y1 )} of a local 
Hamiltonian H on an infinite cylinder of finite width L y , 
where ground state |\f^ yl ) nas anvon charge i inside. This 
approach is valid for any emergent anyon model. Recall 
that if anyon i is Abelian, then di = 1 (as is the case 
for the identity, d\ = 1, which is present in any anyon 
model), whereas if anyon i is non- Abelian, then di > 1. 

For a given ground state |^- yl ), let pi denote the re- 
duced density matrix for half of the infinite cylinder, and 
let {pi, a } denote its eigenvalues. The TEE 7^ for half of 
the infinite cylinder is then obtained from the scaling of 
the entanglement entropy S(pi) = — ^2 a Pi, a l°g(Pi,a) as 
a function of the width L v of the cylinder, which reads 

urn 

S i {L y ) = aL->y i . (DI) 

Following Eq. 40 of Ref. [31] (see also Eq. 18 in 
Ref. [11]), the TEE of half a torus (divided appropriately 
through two cuts in the y direction) in the state \^\ or ) 
with anyon flux of type i in the x direction is equal to 
—2\og(di/D), where the factor 2 comes from counting 
two uncorrelated boundaries. For half an infinite cylinder 
in state |\£^ yl ) we then obtain half of this value, namely 

7i = -log(d i /2?). (D2) 

We can use the numerical estimates of {7^}, and thus 
of {di/D}, for two purposes. First, it follows from the 
definition of D = VEiOi) 2 that T,Mi/D) 2 = 1. There- 
fore, if the sum Y^i{di/D) 2 for our numerical estimates 
is close to one, we obtain a confirmation that the set of 
ground states |\E^ yl ) ^ s complete. If, on the contrary, this 
sum is significantly smaller than 1, this may indicate that 
the set of ground states is not complete. Second, using 
that d\ = 1 and d\ < di, we can use the smallest of the 
estimate {7^} to obtain d\/D, or D. Once we know the 
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total quantum dimension D, we can extract each quan- 
tum dimension di from 7^. 

Notice that a in Eq. IDll is expected to be roughly the 
same for all ground states {|\I/- y1 )}. Indeed, a measures 
the local contributions to the entropy, and all the ground 
states are locally very similar, up to corrections that de- 
cay exponentially with L y . Therefore, for a fixed value 
of L y (e.g. the largest affordable L y in a DMRG calcula- 
tion, which is the one less affected by finite size effects) 
we can also compute di / dj by simply subtracting the en- 



tropies Si and Sj , since Si — Sj 



-{li-lj) = log(di/dj) 



. Then, using that the smallest quantum dimension is 
always d\ = 1, we can compute the rest of quantum 
dimensions di (and therefore also D). This alternative 
approach does not require an extrapolation in L y . How- 
ever, it does not allow us to check that the set of ground 
states {|\I/J y1 )} is complete. 



Appendix E: Computation of the modular matrices 
U and S in a lattice with 71/ 3 rotational symmetry 



Let {|\I/* or )} be an orthogonal basis of (quasi- 
degenerate) ground states of some lattice Hamiltonian 
H on a torus such that H is invariant under a tt/3 ro- 
tation R^/3 (e.g. on a honeycomb, triangular or kagome 
lattice). In this appendix we show how to obtain the U 
and S modular matrices of the emergent anyon model 
from the matrix V of scalar products between the basis 
{|^- or )} and the rotated basis {i^/sl^DL 



Vij = (*T\R 



•tt/3 



\Er tor 

3 



(El) 



Indeed, following Ref. [14], matrix V corresponds to 
the product DWS^D, where D is a diagonal matrix 
with phases Djj = e l€ ^ j (the e l< ^ j corresponds to an arbi- 
trary phase in the definition of the ground state |^ or )), U 

is also diagonal with entries Ujj = e~ l ii c 6j correspond- 
ing to the central charge c of the anyon model and the 
twist 0j of anyon type j, and S has (in general complex) 
entries Sij with Sji — Sij — dj/D > 0. We therefore 
have 



i(Sji) 



* p -i(<f>i-4>j) 



(E2) 



In particular, from 

V u = e-'&e^SuY = e-^Sn, (E3) 

where we used that 0\ = 1 and Su > 0, we can read the 
central charge c (modulus 24). Then, knowing e _z 2T c , 
from 



g 1 24 C 



,(E4) 



for j > 1, where again we used that 0\ = 1 and Sj\ > 0, 
we can read the relative phases e~ t ^ 1 ~ ( ^ j \ with which 
we can compute all e~ % ^ )i ~^ j \ Finally, from 



Vix 



i(Si 



QiS\i, 



(E5) 



for i > 1, where we used that Su > 0, we also learn the 
twist 0i . In this way, we have completely characterize the 
entries Ujj of the modular transformation U. The entries 
Sij of the modular matrix S are then obtained from 



= (SjiY 



(E6) 



where we can use the fact that the exact S matrix fulfils 
Su = Su as a non-trivial check. 

In the above derivation we assumed that we knew 
which state corresponds to identity charge i = 1, and 
placed that state in the first position of rows and columns 
of various matrices. How can we know which ground state 
|^- yl ) corresponds to the identity anyon charge? In a chi- 
ral theory, it might be possible to use the entanglement 
spectrum of the different ground states {|\I/^ y1 )} (or the 
doubled entanglement spectrum of |^- or )) t° identify the 
anyon charge i = 1, for instance as the one containing 
the tower of the identity primary field of the correspond- 
ing CFT, as we did in the main text of this paper for 
the semion model. Alternatively, in some cases the iden- 
tity charge may be the only anyon charge with quantum 
dimension di = 1 (as computed e.g. through the entan- 
glement entropy of |^- yl ), see previous section, or directly 
from the entries of the S matrix). If several states |\£^ or ) 
have the same quantum dimension di = 1 and nothing 
else allows to identify which one corresponds to the iden- 
tity, then the above derivation can be applied repeatedly 
by assuming each time that a different charge with di = 1 
is the identity charge. From the resulting S and U ma- 
trices and the several non-trivial constraints they must 
fulfill, we may then be able to tell whether that assump- 
tion was correct, although we do not know if this will 
always succeed. 



